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ABSTRACT 

VLBI observations at 5 GHz have revealed that supernova 1979C, in the galaxy M100 in the Virgo 
cluster, has shell structure. The shell is approximately circular with the 50% contour deviating from 
a circle by an rms of no more than 8% of the radius. The brightness distribution along the ridge may 
vary. The position of the center of the shell is at R.A.=12 h 22 m 58!66758 and decl.=15° 47' 51'.'9695 
(J2000), with a standard error of 0.8 mas in each of the coordinates. No isolated central component 
is visible above a flux density limit of ^150 /iJy which corresponds to an upper limit of ~15 times the 
corresponding spectral luminosity of the Crab Nebula. The radio lightcurve is clearly falling again 
and the radio spectrum is now flatter than at earlier times. SN 1979C is only the fourth, and oldest 
(optically identified) supernova of which a detailed image could be obtained. 

Subject headings: supernovae: individual (SN 1979C) — galaxies: individual (M100) — radio contin- 
uum: supernovae 



1. INTRODUCTION 

SN 1979C (type ILL) was discovered in the galaxy 
M100 (NGC 4321, di stance: 15.2 Mpc as estimated by 
iFreedman et al.ll2001ft i n the Virgo cluster o n 1979 April 
19 by G. E. Johnson (jMattei et al.l 1 19791 ) and subse- 
quently became one of the most optically and radio lu- 
minous Type II supernovae ever found. On the basis of 
its radio lightcur ve, its explosion da te is estimated to be 
1979 April 4 fsee lWeiler et al.lll986| ). hereafter t = yr. 
VLBI observations started as early as 1982. They re- 
sulted in a determination of the angular expansion rate 
(|Bartellll985t IBartel et all 1 1985ft and deceleration. This 
determination allowed an estimate of the mass-loss to 
wind-velocity ratio of the progenit or, which was found 
to be an order of magnitude lower (jBartel fc Bictcnholz 
[2003J) than that reported from radio lightcurve fitting. 

However, despite a quarter century of VLBI obser- 
vations, the relatively far distance to Virgo has pre- 
vented determining the detailed structure of SN 1979C 
with VLBI. The determination of the expansion rate 
and the deceleration had to be based on size estimates 
made by using a model fit of the supernova bright- 
ness distribution. Clearly, a determination of the su- 
pernova's morphology is important for a more accurate 
determination of the expansion velocity and decelera- 
tion. Such more accurate determination would also al- 
low for a direct, geometric, estimate of the distance to 
SN 1979C and Virgo usin g the expanding shock front 
method (ESM; see, e.g.. IBartel et al.ll2007f) and resolve 
the a mbiguity in an earlier estimate ( Bartel fc Bictcnholz 
2003). It is also important for revealing the mor- 
phological type of the supernova: purely shell-like or 
center filled, perhaps with a bright compact compo- 
nent. Previously, only three other optically-identified 
supernovae could be imaged with an angular resolution 
high enough to unambiguously dete rmine the morpho- 
logical type. All three, SN 1986J (jBartel et all Il991t 
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Bietenholz e t alj |2004ft SN 1987A dGaensler et al.lll997l 
2007ft . and SN 1993 J (iMarcaide et ; alJll995t IBartel et al l 
19951 l200d : iBietenholz et all 12003ft have shell structure, 



with SN 1986 J having a central component in addition 
to the shell. SN 1979C is the oldest optically-identified 
supernova that can still be monitored with VLBI. As 
such, it provides a bridge between our understanding of 
supernovae (interacting with the circumstellar medium 
left over from the progenitor) and supernova remnants 
(interacting with the interstellar medium). 

Only a few supernova remnants which are less than 
a ce ntury old are known, all extragalac tic, e.g., in 
M82 ([McDonald et al.ll2001ft and in Arp220 ([Parra et all 
12007ft . The youngest Galactic remnants are several 
centuries old, e.g., the Crab Nebula, which shows 
a pulsar wind nebula, Cas A, which shows a shell, 
and G21.5 — 0.9, which shows both a pulsar wind neb- 
ula (e.g., [Bietenholz fc Bartell l2008bft and likely also 
a sh ell ([Matheson fc S afi-Harb! [20051 : iBocchino et al.l 
12005ft . Since all radio supernovae originate from the ex- 
plosions of massive stars, they are expected to leave a 
neutron star or a black hole. Either a young neutron 
star or a black hole may be accompanied by a bright, 
compact radio source near the center of the expanding 
shell. Young neutron stars are likely to be associated 
with radio-bright pulsar wind nebulae, and black holes 
may have radio bright jets due to an accretion-disk sys- 
tem. 

Indeed, pulsar wind nebulae are found around a num- 
ber of young pulsars, although the youngest known pul- 
sars have ages near 1000 yr. A central radio source pos- 
sibly associated with a neu tron star or black hole was re- 
cently found in SN 1986 J ([Bietenholz et alj|2004ft . Such 
a central radio source might also be seen near the cen- 
ters of other older supernovae when clearly resolved im- 
ages are obtained. Here we report on such an image 
of SN 1979C, discuss its morphology, radio lightcurve, 
and spectral evolution, and search for signs of a central 
compact remnant of the explosion. The expansion rate, 
deceleration, and distance determination will be reported 
elsewhere. 
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2. OBSERVATIONS AND DATA ANALYSIS 

The VLBI observations were carried out on 2005 
February 25 with a global array of 21 antennas. The 
array consisted of the Very Long Baseline Array (VLBA; 
each of 25 m diameter), phased Very Large Array (VLA; 
130 m equivalent diameter), and the Robert C. Byrd 
telescope (GBT; ^105 m diameter), all of the National 
Radio Astronomy Observatory, NRAO, in addition to 
the Arecibo Radio Telescope (305 m diameter) and eight 
telescopes of the European VLBI Network, namely the 
Effelsberg Radio Telescope, Germany (100 m diameter), 
the Westerbork Array The Netherlands (92 m equiva- 
lent diameter), and the telescopes at Medicina and Noto, 
Italy (32 m diameter), at Torun, Poland (32 m diameter), 
Jodrell Bank, Great Britain (26 m diameter), and Hartc- 
beesthoek, South Africa (26 m diameter). We observed 
for a total of ~15 h at a frequency of 5 GHz with a band- 
width of 32 MHz and recorded both senses of circular po- 
larization with a total recording bit rate of 256 Mbits s -1 . 
The data were correlated with the NRAO VLBA proces- 
sor in Socorro and further analyzed with NRAO's As- 
tronomical Image Processing System (AIPS). The ini- 
tial flux density calibration was done through measure- 
ments of the system temperature at each telescope, and 
then improved through self-calibration of the reference 
sources. 

We phase-referenced our VLBI observations to 
J1215+1654, which is 2.2° away on the sky, and used, 
through most of the run, a repeating cycle of 250 s, point- 
ing for 170 s to SN 1979C and for 80 s to the reference 
source. Finally, we also used the VLA, which was in the 
B configuration, to determine the total flux density of 
SN 1979C at 1.43, 4.99, and 8.43 GHz. In order to ob- 
tain the best possible image, we phase self-calibrated the 
data after phase-referencing. 

3. RESULTS 
3.1. The Brightness Distribution 

In FigureCQwe show the VLBI image of SN 1979C jux- 
taposed with an earlier VLA image of the host galaxy 
M100. The supernova is located in a southern spiral arm 
of the galaxy. The VLBI image shows that the supernova 
has shell structure with a brightness minimum located 
only slightly to the west of the center. The outer con- 
tour at half of the maximum brightness is approximately 
circular with rms deviations from circularity of 8% of the 
radius. We think that this estimate is somewhat conser- 
vative since it is based on an image convolved with an 
elliptical beam, which tends to elongate the structure 
along the long axis of the beam, which is oriented north- 
south. The 50% contour is indeed slightly elongated in 
the north-south direction. An ellipse fit to the 50% con- 
tour has a ratio between the major and minor axes of 
1.10 and a p.a. of 173°. 

The maximum brightness is only 187 /iJy beam -1 
and the rms noise brightness in the background is 
11 /xJy beam -1 . How much could possible deconvolu- 
tion errors contribute to the uncertainty in the bright- 
ness distribution of the supernova? We investigated the 
significance of the appearance of the shell structure and 
of features in the brightness distribution by using simula- 
tions. We used a disk model as a test source, converted 
it to visibility data for our u-v coverage, added noise, 



and then imaged and deconvolved the simulated data in 
the same way as the observed data. In this way we pro- 
duced 36 images of our model source with different noise 
realizations but with the rms background noise bright- 
ness being always close to our observed background rms 
of 11 /xjy beam" 1 . The average of the 36 images resem- 
bled the original disk model, indicating that there are no 
systematic distortions or biases. The individual images, 
however, varied in their appearance with the rms of the 
brightness variations on the image being twice as large 
as the rms of the background noise brightness. None of 
these images had a closed ridge with a clear minimum 
in or near the center. We therefore exclude the possibil- 
ity that the true brightness distribution is uniform across 
the source or has a filled center. Instead, with high prob- 
ability, it is indeed a shell. 

Our simulations also showed that the northern and 
southern low-brightness extensions are likely due to cor- 
rugations in the image plane caused by our incomplete 
u-v coverage, and are therefore probably not signifi- 
cant. The brightness varies along the ridge between 
187 /j,3y beam -1 and 112 /iJy beam -1 , representing a 
range of only ±1.7x the estimated on-source brightness 
uncertainty of ~ 22 /zJy beam -1 . It is not clear whether 
the variation is significant. 

The brightness at the local minimum slightly to the 
west of the center of the shell is 80 /iJy beam -1 . 
The brightness increases toward the east and is ~ 
130 [iJy beam -1 at the center, although no isolated cen- 
tral component is visible. We determine an upper limit of 
the flux density of any compact component in the central 
area of the shell of ~ 150 fiJy beam -1 . 

3.2. The Position 

Phase-referencing allowed us to determine the posi- 
tion of SN 1979C relative to that of our calibrator, 
J1215+1654. The calibrator's posit ion was taken from 
the VLBA Calibrator Survey fVCSl. lBeaslev et al.ll200l 
as R.A.=12 h 15 m 03!979123 and decl.=16° 54' 37/95733 
(J2000). The standard error 3 is 0.5 and 0.6 mas in 
RA and deck, respectively. In our VLBI image of 
SN 1979C in Figure HJ we place the origin of the co- 
ordinate system at the estimated center of the outer 
50% contour. This coordinate origin is at the position 
of R.A.=12 h 22 m 58!66758 and decl.=15° 47' 51'.'9695 
(J2000). We estimate that our standard error of the po- 
sition of SN 1979C relative to J12 15+1654 is 0.5 ma s 
in each coordinate (compare with iPradel et al.l [2006) . 
Adding in quadrature the standard position error of 
J1215+1654's to that of SN 1979C, we estimate a stan- 
dard error of the position of the center of the image in 
Figure [T] to be 0.8 mas in each coordinate. 

As a check, we also determined the position of 
SN 1979C relative to J1215+1654 by using our earlier 
1.7 GHz data from 2001. At this frequency, SN 1979C 
was barely resolved, and contributions to the position er- 
ror due to unmodeled effects of the ionosphere expected 
to be much larger. We estimate a standard error of 2 mas 
in each coordinate. Our estimated relative position of 
SN 1979C at 1.7 GHz in 2001 was within 1 mas of the 
more accurate position given above, well within the un- 
certainties. 

3 IBeaslev etail (200% call this error "inflated formal error." 
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3.3. VLA Total Flux Densities 

We determined the total flux density, S v , of SN 1979C 
from the interferometric observations with the VLA 
made during the time of our VLBI observations. We 
list these flux densities for frequencies, u, of 8.43, 4.99, 
and 1.43 GHz in Table [T] In addition, we give the flux 
densities for v — 8.46 and 1.67 GHz measured in 2001, 
for comparison. Then we convert the flux densities to 
v = 4.89 and 1.47 GHz and plot them in Figure 2, to- 
gether with the much larger s et of flux density determ i- 
nations at these freq uenci es oflWeiler et alj (119861I1991D , 
iMontes et all (|2000D . and lBartel fc Bietenholzl (|2003| ). 

The new measurements of the flux density show that 
the supernova is again fading, reaching at 1.5 GHz its 
lowest level ever. The bumpy pattern seen previously is 
extended downward. 

3.4. Radio Spectra 

IMontes et al.l (|2000D found that the optically-thin spec- 
trum had a spectral index, a, of —0.75 with a deviation 
range of -0.76 to -0.63 for t < 12 yr (i.e., before 1991; 
a being obtained from the fit of an empirical model to 
the multi- frequency radio lightcurves) , and a ~ —0.70 for 
t > 12 yr. At their latest times, t = 17.7 and 18.8 yr, they 
obtained a = —0.63 ± 0.03. Our latest measurements 
give still flatter spectra (see Table [T}. The difference 
of our spectral indices to that of —0.63 may be signifi- 
cant. In Figure [3] we plot our flux density measurements 
and the corresponding spectra for our two latest epochs. 
For comparison we show the f lux density measurem ents 
at t = 17.7 and 18.8 yr from IMontes et ail (|2000l ) and 
the corresponding average spectrum. Clearly, there is 
evidence that the spectrum of SN 1979C has flattened 
further since t ~17 yr. 

4. DISCUSSION 

SN 1979C is the fourth (optically identified) extra- 
galactic supernova for which an image with sufficient an- 
gular resolution could be made to determine the nature 
of the brightness distribution. The other three, as men- 
tioned in the introduction, are SN 1986J, SN 1987A, and 
SN 1993J. Our VLBI image shows that SN 1979C has 
a morphology similar to the other three, namely a fairly 
to highly round shell, with the brightness possibly be- 
ing modulated along the ridge. The possible brightness 
modulation likely arises from density irregularities either 
in the circumstellar medium or in the ejecta. 

The flux density of SN 1979C is now clearly again 
on the downward trend. Earlier, it was suggested that 
the radio lightcurv e exhibits a quasi-periodic pattern 
(|Weiler et al.l [1992) which was interpreted as possibly 
being due to the dynamical modulation of the wind 
from the red supergiant proge nitor by a binary compan- 
ion (|Schwarz fc Prindd fl996h . With more data up to 
t ~ 19 yr, it was noticed that any periodic variation of 
the radio lightcurve was not extended beyond t ~ 10 yr. 
Instead, it was suggested that, with the flattening of the 
radio lightcurve, the su pernova had apparen tly entered a 
new stage of evolution (|Montes et al.ll20doh . 

With our new data up to t ~ 26 yr, which show that 
the lightcurve is again falling after t ~ 16 yr, a new 
look at the characteristics of the radio lightcurve is war- 
ranted. In particular, we think that the evidence for 



a periodicity in the radio lightcurve or a new stage of 
evolution now seems weaker. If at all, then the radio 
lightcurve appears to be periodic on a logarithmic time 
scale. However, a closer inspection leads us to think that, 
rather than being periodic on a logarithmic time scale, 
the radio lightcurve is modulated on a whole range of 
time scales. For instance, there appear to be variations 
on time scales of less than a year. There are also vari- 
ations on time scales of several years with e.g., a 5-yr 
bump in the radio lightcurve at t ~ 7 yr and a bump 
lasting at least 15 yr at the end of the observed radio 
lightcurve. These variations are likely caused by varia- 
tions in the density o f the ejecta (for discussion s in the 
case of SN 1993J, see iMioduszewski et al]|2001ft and/or 
of the circumstellar medium. In case of the variations 
in the latter, the radio lightcurve can be thought of as 
a "time machine" which records in reverse the mass-loss 
to wind-velocity ratio of the progenitor for thousands of 
years before the star died. 

The spectrum appears to have flattened somewhat 
with the spectral index being above the dev i ation range 
of -0.76 to -0.63 given by IMontes et all (|2000h . A 
similar behavior is also seen for SN 1993 J where the 
spectrum flattens while the source dimming accelerates 4 
(|Bartel et al.l 120021 ; IPerez-Torres et al.ll2002l) . 

No clearly isolated component was found in the cen- 
tral region of the projected shell of SN 1979C. How- 
ever, the brightest portion of the shell appears to ex- 
tend toward the center in projection, and could overlap 
with any compact component that may be present in 
the interior. An isolated central component was in fact 
found for SN 1986J, with a spectrum peaking at ^20 GHz 
while the emission from the remaining parts of the su- 
pernova ha s an optically thin spectrum from at least 1.7 
to 22 GHz (jBietenholz fc Barte]|[2008al : iBietenholz et alj 
l2004f ). This isolated central component may be emission 
associated with the compact remnant of the explosion. 

It is intriguing that the spectrum of SN 1979C shows 
evidence of flattening with time. Such flattening would 
be expected for a supernova with emission from a com- 
pact remnant that is starting to appear in the center 
of a shell which is expanding and becoming increas- 
ingly transparent at progressively higher frequencies. 
However, no central component can clearly be seen in 
SN 1979C, so the observed spectral flattening cannot be 
ascribed to such a component. 

It is probable, therefore, that the spectral flattening 
is intrinsic to the shell's synchrotron emission. The flat- 
tening must then be due to a change in the population 
of accelerated particles which produce the synchrotron 
emission. Calculations have shown that physical condi- 
tions in supernova s hocks can effect the spect ral index of 
radio emission (e.g., lEllisonlt~aT][2007L l200l . 

Our new upper limit on the flux density of any cen- 
tral compact source in SN 1979C of 150 /xJy puts 
strong constrai n ts on the luminosity of such a source. 
iBandiera et alj (|1984l ) argued that young pulsar wind 
nebulae would have spectral luminosities of 10 to 1000 
times that of the Crab Nebula. The spectral lumi- 

4 Please note that although Wcilcr et al. (2007) claim that the 
spectral index of SN 1993J remains constant over time, an in- 
crease of the spect ral in dex consistent with the m easurements of 
IBartel et all J2002T) and IPerez-Torres et alj J2002T) can be clearly 
seen in their data with the smallest errors and/or scatter. 
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TABLE 1 

Total flux densities and spectral indices of SN 1979C 



Date 


Age (yr) a 


1.43 GHz 


Flux densities (mjy) b 
1.67 GHz 4.99 GHz 


8.43 GHz 




2001 Feb 23 
2005 Feb 25 


21.89 
25.90 


3.19±0.22 


4.28±0.23 

1.68±0.09 


2.3±0.4 dc 
1.37±0.10 


-0.38 ±0.15 
-0.49 ± 0.09 



a Time since assumed explosion date of 1979 April 4 (1979. 26). b Total flux density, S v , measured with the VLA. The uncertainties 
are approximate standard errors and were computed by adding in quadrature the statistical standard error and a systematic 
calibration error of > 5%. c The spectral index, a, determined from the flux density measurements, S v , at frequency, v: S v cx v a . 
For our last epoch a weighted fit was use d. The errors are one s tanda rd deviation, determined by varying the flux density values 
by their standard errors. d Data are from lBartel fc Bietenholz! (2003).° The observing frequency was 8.46 GHz. 



nosity of the central component of SN 1986J is ~200 
times th at of the Crab Nebula bet ween 14 GHz and 
43 GHz (|Bietenholz fc Bartelll2008al) . In contrast, the 
upper limit on the 5 GHz spectral luminosity of any 
central component in SN 1979C is only 15 times the 
corresponding spectral luminosity of the Crab Nebula, 
more than an order of magnitude fainter than both 
that in SN 1986J an d our previous limit for SN 1979C 
dBartel fc Bietenholz! 120051) , and at the low end of 
iBandiera et alj ( 1984f )'s estimate. Any pulsar wind neb- 
ula that may exist inside the shell can probably not yet 
be detected because the shell is likely still opaque at ra- 
dio frequencies. Since the supernova shell becomes less 
dense and more transparent as it expands, any central 
component will become more easily detectable in the fu- 
ture. Searching for emission from the compact central 



remnant of the explosion therefore remains an intriguing 
prospect. 
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Fig. 1. — Left side: A radio image of the spiral galaxy M100 (NGC 4321) in the Virgo cluster, with SN 1979C in it, made with the 
VLA in the C array at 1.7 GHz on 1996 March 15. The contour levels are at 1.3, 2.5, 4, 8, 16, 32, and 64% of the peak brightness of 
15.8 mjy beam" 1 . The ful l-width at half- maximum (FWHM) restoring beam is plotted in the lower left corner. The image is taken from 
IBartel fc Bietenholzl J2003T) . Right side: the VLBI image of SN 1979C, made at 5.0 GHz on 2005 February 25. The contour levels are at 
— 17, 17, 30, 40, 50 (thick white contours), 60, 70, and 90% of the peak brightness of 187 /xjy beam -1 . The grey scale, in /xjy beam -1 , 
is given at the top. The rms brightness of the background noise is 11 (j,Jy beam -1 . The FWHM of the restoring beam is 2.0 mas X 
1.1 mas oriented at a p. a. of —10°. North is up and east to the left. The origin of the coordinate system is at R.A.=12 h 22 m 58?66759 
and decl.=15° 47' 51'.'9695 (J2000) with a standard error of 0.8 mas in each coordinate. It is at the approximate center of the outer 50% 
contour (see text). 



6 



Bartel and Bietenholz 



20 



w 
£ 
0) 

X 



10 



4 



a, a 1.47 GHz 
o,m 4.89 GHz 



4 6 8 10 20 

Time since explosion (yr) 



Fig, 2, — The r adio lightcurves at 1.47 and 4.89 GHz. Open symbols rep r esent values taken from Wcilcr ct al. (1986, 1991]) and 
IMontes ct al. (2003). Filled symbols represent values from Bartel & Bietenholz (2003) and this paper. The observed valu es were taken 
at a number of slightly different frequencies. If necessary, they were converted to the above frequencies with a = —0.75 (Mo ntes et al.l 
2000), except for our values at t = 22 and 26 yr, for which we used the spectral indices from Table [T] The 4.89-GHz value at t = 22 yr is 
interpolated between the 1.67 and 8.43 GHz values from Tablc[T] 
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Fig. 3. — The flux densities from Table [T] for the e pochs given in the lo wer left. For comparison and to show evidence for the flattening 
of the spectrum with time, the flux densities given by Montes et al. (2000) for the time range from t = 17.68 to 18.82 yr after the explosion 
are also plotted. The straight lines indicate fits to the corresponding data. 



